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Warm dark matter (WDM) may resolve the possible conflict between observed galaxy halos and the halos 
produced in cold dark matter (CDM) simulations. Here we present an extension of MSSM to include WDM 
by adding a gauge singlet fermion, % w hh a portal-like coupling to the MSSM Higgs doublets. This model 
has the property that the dark matter is necessarily warm. In the case where M% is mainly due to electroweak 
symmetry breaking, the % mass is completely determined by its relic density and the reheating temperature, 
T R . For 10 2 GeV < T R < 10 5 GeV, the range allowed by % production via thermal Higgs annihilation, the % 
mass is in the range 0.3-4 keV, precisely the range required for WDM. The primordial phase-space density, Q, 
can directly account for that observed in dwarf spheroidal galaxies, Q ~ 5 x 10 6 (eV/cm 3 )/(km/s) 3 , when the 
reheating temperature is in the range Tr ~ 10 — 100 TeV, in which case M% ~ 0.45 keV. The free-streaming 
length is in the range 0.3-4 Mpc, which can be small enough to alleviate the problems of overproduction of 
galaxy substructure and low angular momentum of CDM simulations. 

PACS numbers: 12.60Jv, 98.80.Cq, 95.35.+d 



Introduction: Cold dark matter (CDM) with a cosmologi- 
cal constant (A CDM) is remarkably successful in explaining 
the large scale structure of the observed universe. Numerical 
simulations based on the ACDM model predict cusped central 
densities |Q]]. The observational situation is less clear. It has 
been suggested that observed galaxy halos have cores [2, 3], 
implying that CDM has a cusp problem. However, the obser- 
vational evidence has also been interpreted to support cusped 
halos [4]. Here we will consider the former possibility. One 
way to solve the cusp problem, should it exist, is to have dark 
matter with a sufficiently large velocity dispersion, known as 
warm dark matter (WDM) [5]. In addition, CDM simulations 
result in overproduction of galactic substructure [6] and low 
angular momentum [7], which may be alleviated by reducing 
power at small scales via free-streaming of WDM. 

Several candidates have been suggested to account for 
WDM, notably sterile neutrinos in a minimal extension of the 
Standard Model (MSM) Hi], and superWIMPs d. These 
models may also account for other phenomena, such as pulsar 
velocities and baryogenesis in the case of sterile neutrinos, 
and can have distinctive collider phenomenology, as in the 
case of a gravitino superWIMP, which requires large masses 
> 500 GeV for the NLSP of the MSSM. However, the mass 
of the dark matter candidate in these models is not fixed by 
the model, in which case the dark matter is not necessarily 
warm. The unique feature of the model we will present here 
is that the dark matter candidate is fixed by its relic density to 
be necessarily warm. 

In the case of the MSSM, the only possible WDM candidate 
is a stable gravitino [11]. (Sterile neutrinos could also play 
this role in the MSSM extended to include neutrino masses.) 
However, this is not compatible with the MSSM in the case of 
gravity-mediated SUSY breaking as the gravitino is too heavy. 
The form of SUSY breaking may be determined at the LHC 
from the pattern of SUSY particle masses. Here we present a 



new WDM candidate, the Zi-singlino, in a portal-like exten- 
sion of the minimal supersymmetric (SUSY) standard model 
(MSSM) (HE! in which the stability of the WDM particle 
is ensured by a Z2-parity. The Z2-singlino was previously in- 
troduced to provide a stable dark matter candidate in R-parity 
violating SUSY models 11211 . Here we will show that, in a 
different region of its parameter space, the model provides a 
dark matter particle which is necessarily warm and with the 
right properties to explain non-singular galaxy halos and to 
alleviate the galaxy substructure and low angular momentum 
problems of CDM models. 

The Model: We extend the MSSM by adding a chiral su- 
perfield % and messenger field S of mass M$. We also impose 
an additional Zz symmetry under which % is odd, while all 
other fields are even. The effective superpotential after inte- 
grating out S is given by lfl2ll 

W = W M ssm-\ — 1 — . (1) 

Ms I 

Since % is odd under Z2, its lightest component cannot decay 
to any of the MSSM fields. Therefore the fermionic compo- 
nent of %, the Z2-singlino, % is a good candidate for DM . In 
Eqn. ([T} we have included a SUSY mass M% Q for %. Since in 
general the % mass must be small relative to the weak scale 
in order to account for WDM, a particularly interesting case 
is where the % mass is entirely due to the Higgs expectation 
values, with M^ o zero or negligibly small. The small % mass 
can then be understood in terms of a large value for M$ com- 
pared with the Higgs expectation values. The absence of a X 
mass term in Eqn.lfl]) is guaranteed if there is an unbroken R- 
symmetry which allows the /jH u H £ / term of the MSSM, since 
in this case the R-charge of % must be zero. 

Z2-singlino dark matter is interesting as a SUSY implemen- 
tation of gauge singlet dark matter. Gauge singlet scalar dark 
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matter interacting via the Higgs portal |[l4ll was first proposed 
in [15] for the case of complex scalars and in for real 
scalars, and was further discussed in OHdEoElll. Cou- 
plings to hidden sector particles are currently of considerable 
interest JH |H |H |H |2^, HI, IH |H|. The simplicity of the 
terms in Eqn.([T]) should allow them to easily form part of a 
hidden sector dark matter model. 

In the following we will calculate the relic abundance of % 
as a function of its mass and reheating temperature, Tr. We 
will show that in the case where the % mass is mostly due 
to the Higgs expectation value, % is necessarily warm with 
a mass in the keV range when it accounts for the observed 
dark matter density. For reasonable values of Tr the primor- 
dial phase-space density can then account for that observed in 
dwarf spheroidal galaxies (dSphs), while the free-streaming 
length can damp the density perturbation on small scales and 
so reduce galaxy substructure formation and angular momen- 
tum loss. 

Relic Abundance of %: Production of % will occur mainly 
through thermal Higgs annihilations 2 . We will see that most 
of the x are produced at temperatures close to the reheating 
temperature, Tr. In this case we can consider the Higgs ex- 
pectation values to be zero and calculate with the weak eigen- 
state Higgs doublets. % production via thermal Higgs boson 
pair annihilation occurs due to the contact interaction in the 
Lagrangian ,fxxH u H c i /Ms- The total rate of % production per 
Higgs pair annihilation is then 



dtij 
dt 



+ 3Hn 1L = T 1 rn H ; Tyr = 8n H a H v re i 



(2) 



In this Tin = 2AT 3 /% 2 is the number density of a complex 
Higgs scalar, Oh — {f/Ms) 2 /16n is the cross-section for 
h° u h° d — > xx (with the same for h+hj — > XX) an d v rel — 2 is the 
relative Higgs velocity. There is an overall factor of 2 in 
since each Higgs pair annihilation produces two % particles. In 
addition, there will be % production via Higgsino-Higgs boson 
annihilation to XX°- This will increase the % production rate 
by approximately a factor of 2, which we have included in Fj. 
In terms of T we obtain 



2.4M 



7tVsin 2 2p ' 



(3) 
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where tanP = v„/v c / with v = yv^ + v 2 , 
is the contribution to the x mass from electroweak symme- 
try breaking, Mj s y = /v 2 sin2p/M,$ . This relates the sym- 
metry breaking contribution to the % mass to the strength of 
the % interaction with the MSSM Higgs, f/M s . As a re- 
sult, the % relic density fixes the x mass in the case where 
it is mostly due to symmetry breaking. The total x mass is 
then A% = +M^ sb . Since Yj <x T 3 while H oc T 2 during 



radiation-domination, the production of % will occur mostly at 
the highest temperature during radiation-domination, which 
is the reheating temperature, Tr. (The temperature can be 
higher during the inflaton-dominated era before reheating, but 
since H T 4 during this era, its contribution to x produc- 
tion is small compared with production at reheating.) We 
define the thermalization temperature 7J/j by the condition 
r f = H = k T T 2 /M Ph where k T = (4n 3 g(T)/45)^ 2 and g(T) 
is the number of relativistic degrees of freedom. Therefore 



T t h = 



/t r 7i;Vsin 2 2p 
2AM 2 sb M PI 



Assuming g(T) is constant, Eqn. (O can be written as 

1 n H 



dT \ T 3 J T th T 3 ' 



(4) 



(5) 



where we have used the relation F%/HT = 1/7^, which is 
generally true during radiation-domination. Integrating this 
from Tr to T, and noting that nn/T 3 = 2.4/ft 2 is a constant, 
we obtain 



T 3 



2A(T R -T) 
n 2 T th 



(6) 



Clearly most of the production of the comoving % density 
occurs at T w Tr. Including an entropy dilution factor for 
the change of g(T) from Tr to the present CMB temperature, 
Ty = 2.4 x 10~ 13 GeV, we find for the present x number den- 
sity, 



n x(Ty) 



2AT 3 g( Ty ) T R 
n 2 g(T R )T th 



(7) 



Using Eqn.© to eliminate 7J/,, the relic abundance of x from 
Higgs annihilations is then given by 



4(1.2) 2 (M lsh \ 2 g(Ty) T R T 3 M 3 M PI 



k 5 \ Mj ) g{T R ) ^y4 s in 2 2p p c - 



(8) 



where p c = 8.1 x 10 h GeV is the critical density. There- 
fore 



My = 1.92 



M %sb 



'_ \ 2 / 3 



•lTeVX 1 / 3 . 2/3 



0.113 



2pkeV 



2\ 1/3 



(9) 



where we have used g(Ty) — 2 and g(T p ) = 228.75, corre- 
sponding to the MSSM degrees of freedom, and we have 
expressed £1^ relative the observed dark matter abundance, 
%/i 2 = 0. 1 1 3 1 ± 0.0034 JH - 

From Eqn. (0 we see that in the case where the x mass is 
due to symmetry breaking, Mj is automatically close to a keV. 
As a result, the dark matter in this model is necessarily warm. 
The relationship between M% and Tr for different values 

will therefore neglect the contribution of Higgs decays in the following. of tanp is shown in Fig Q] for the case Mj s h = Mj. From 



2 X can also be produced by decay of thermal Higgs particles. However, 
we find that this is negligible compared with Higgs annihilation except for 
300 GeV, where both processes give a similar contribution to nj. We 
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FIG. 1 : Contours of %ft 2 = 0.1131 ± 0.0034 are shown in the plane 
of Mjr versus Tr for different values of tan|3. The % mass is set to its 



value from symmetry breaking, Mj = M- 



of view, the two models are not dissimiliar in their need for 
tuning. Such tuning is a generic problem for any dark mat- 
ter model which does not directly address the baryon-to-dark 
matter ratio. 

The Phase-Space Density of %: Dark matter with a finite 
phase-space density may be able to explain the finite density 
of galaxy cores [32]. The coarse-grained phase-space density 
is defined by Q = p%/o^., where the 1-D velocity dispersion 

is = ^(1/3)(P 2 /M 2 ) [3Z] . If Q is finite then there is a 

limit on how large p can be for a given velocity dispersion, 
which prevents the formation of singular galaxy cores. Q can 
be expressed in terms of the distribution of % produced by 
thermal Higgs annihilation i33ll 



Q 



3 3/2 ^4Px = 3 3/2 Mf% p c J 

(pi?' 2 n 



3/2 
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this we find that the range of % mass is tightly constrained. 
Tr ^ 10 6 GeV is necessary in order to avoid thermal gravitino 
overproduction (for the case with hadronic gravitino decay 
modes) ll3~ill while Tr > 10 2 GeV in order to have the ther- 
mal Higgs necessary to produce the % density. In the case 
where the % mass is due to spontaneous symmetry breaking, 



M-- 



M%, the % mass is in the range 



0.19 sin 2/3 2p keV<A%<4.1 sin 2/3 2p keV 



(10) 



Moreover, Tr < T t h is necessary in order to keep % out of ther- 
mal equilibrium, which is assumed in Eqn.©. This gives 
My ^,0.3 keV. Thus for reasonable values of sin2p we ex- 
pect 0.3 keV < Mj < 4 keV. The corresponding range of re- 
heat temperature is 10 2 GeV < T R < 10 s GeV. 

In this we have used the observed abundance of dark mat- 
ter as an input, which then determines the mass of the dark 
matter particle to be of the order of a keV. In doing so, we 
have implicitly tuned the mass of the dark matter particle for 
a given Tr. To put this tuning in perspective, we can compare 
the present model with the most popular dark matter scenario, 
that of thermal relic weakly interacting dark matter. In this 
case, the dark matter density is broadly of the correct order of 
magnitude when the dark matter particle mass and its inter- 
action strength are determined by the weak scale. However, 
although phenomenologically encouraging, theoretically this 
amounts to an implicit tuning of the physics of thermal freeze- 
out against the completely unrelated physics of baryogene- 
sis in order to obtain abundances of baryons and dark matter 
which are within a factor of 6. In other words, the interaction 
strength, which is determined by its mass scale, must be tuned 
in order to obtain the correct dark matter abundance relative 
to the baryon density. Similarly, in our model we are tuning 
the interaction strength, which is now determined by the dark 
matter particle mass, to obtain the correct dark matter abun- 
dance relative to the baryon density. In the thermal relic case 
the output is an interaction mass scale of the order of the weak 
scale. In our case the output is a mass for the dark matter par- 
ticle of the order of a keV. Therefore, from a theoretical point 



where T c { — (g(T y )/g(TR)) l ^Ty for a decoupled density cre- 
ated at T R and J = £y 2 f{y)dy/ S^y A f{y)dy. Here f(y) is the 
distribution function for production of % from thermal Higgs 
annihilation, where y = \P%\/Td is the comoving momentum. 
Therefore 



Q = 5.6x 10 7 — 
* I 0.113 



1 keV 



J 

oT 



3/2 



eV/cm 3 



(km/s) 3 



(12) 

The distribution function due to thermal Higgs annihilation 
is not known at present. However, J may be estimated by 
assuming that the momentum of the % at Tr is equal to the 
mean momentum of the thermal Higgs particles, p « 3T, in 
which case f(y) — 8(y — 3). This gives 7=1 jy 1 = 1/3 2 = 
0.11. This is in good agreement with the value obtained in the 
case of thermal Higgs decay, J — 0. 12 ll34il . 

Q is conserved for an adiabatically contracting collision- 
less gravitational system j32tl . More generally, Q can only 
decrease from its primordial value during relaxation in col- 
lisionless systems. Simulations of structure formation have 
shown that the phase-space density can decrease by a factor 

10 2 to 10 3 iH. 

If the distribution of dark matter in dwarf spheroidal galax- 
ies (dSphs) is cored, then a synthesis of recent photometric 
and kinematic dSph data [3], which indicates a mean density 
~ 5GeV/cm 3 and central velocity dispersion ~ lOkm/s, im- 
plies that HI 



QdSph «5xlO*^f 
(km/s) J 



(13) 



Eqn. (fT2l then provides an upper limit on My_ in the case with 
an explicit SUSY % mass. If we assume a dynamical sup- 
pression of Q by at most a factor 10 3 , then the largest value 
of Q consistent with observations is ps 10 10 eV/cm 3 /(km/s) 3 . 
From Eqn.(Q~2]l the largest mass compatible with this is M% w 
6 keV. This is not much larger than the mass range 0.3-4 keV 
implied by the relic density in the case with Mj w s i,. Since 
there is no reason to expect a SUSY mass for % in the keV 



4 



range, it is more natural to assume that Mj Q is zero or neg- 
ligible, with the small % mass then generated by electroweak 
symmetry breaking together with Ms 3> v. 

Using Eqn. (0 to eliminate Mj from Eqn.([T2b, the phase 
space density of % as a function of Tr is given by 



Q=4.0x 10 8 sin 2 2(3' x 



M lsb 



0.113 



1 TeV 

Tr 



J 



3/2 



eV /cm 3 



0.1 J (km/s) 3 



(14) 



Comparing with Eqn. ( fT3l . we see that the value of Q can 
be equal to or larger than 2dSph f° r reasonable values of Tr. 
From Eqn. ( fT4b the required reheating temperature is 



r fl «80 sin 2 2(3 



6dSph 



My 



M- 



1 & 



TeV 



(15) 



Thus in the case where the % mass is mostly from the Higgs 
expectation value and the phase-space density of dSph cor- 
responds to the primordial phase-space density without sup- 
pression, 10 TeV "zTr"Z 100 TeV is typically required, de- 
pending on sin2p. From Eqn.© the corresponding % mass 
is M^ ~ 0.45 keV. If the primordial phase space density 
is suppressed during formation of dSphs, then larger Q and 
lower Tr are required. In general, for sin 2(3 > 0.1 and 
10 2 GeV < T R < 10 s GeV, Q due to % dark matter is in the 
range 10 5 - 10 10 eV/cm 3 /(km/s) 3 . 

Free-streaming length of %: The free-streaming length, 
A,f s , below which primordial perturbations are suppressed, is 
roughly equal to the horizon when the % particles become non- 
relativistic. In general, for a distribution of relativistic decou- 
pled particles, Xf s is given by ll36ll 



Xf S ks 1.2 Mpc 



1 keV 

My 



10-75\ 1/3 ( (p/T) 
s{Tr) ) \ 3.15 



(16) 



where (p/T) is the mean momentum over T of the initial rel- 
ativistic % distribution. Since the % are produced by annihi- 
lation of thermal Higgs at T Tr, we expect that their mean 
momentum will be approximately equal that of the thermal 
Higgs at T R , such that (p/T) w 3. Thus with g(T R ) = 228.75 
and 0.3 keV < Mj < 4 keV, X fs is in the range 0.3 - 4 Mpc. 
(Smaller values are possible if the mean momentum of the % 
distribution is less than thermal.) The lower end of this range 
is comparable with the scale of galaxies and so may alleviate 
the problems of overproduction of substructure and low angu- 
lar momentum observed in CDM simulations 

Lyman-oc constraints: Observation of Lyman-a absorp- 
tion spectra constrains the matter power spectrum on small- 
scales and so provides a lower bound on the WDM par- 
ticle mass given its momentum distribution l37tl . Lower 
bounds on the mass of sterile neutrino WDM were obtained 
in JM El S3). In 0], lower bounds were obtained for 
WDM momemtum distribution functions which are general- 
izations of the Fermi-Dirac distribution. In the case of thermal 



relics which decoupled while relativistic, a lower bound of 1 .7 
keV (95% c.l.) was obtained, while for non-resonantly pro- 
duced sterile neutrinos (less-than-equilibrium density but with 
a thermal momentum distribution) the corresponding lower 
bound was9.5keVll4lll. 



Since the distribution function from thermal Higgs anni- 
hilation is not expected to be of Fermi-Dirac form, we can- 
not directly apply the results of existing Lyman-a analyses 3 . 
However, as the % density from thermal Higgs annihilation is 
less than the thermal equilibrium density but the mean % mo- 
mentum is of the order of the thermal Higgs momentum at 
T « niH, we expect the Lyman-a lower bound on the % to lie 
between the lower bounds of |4lll . which may allow a window 
below the 4 keV upper bound from the % relic density. 

In addition, Lyman-a observations can easily be consistent 
with WDM in the case where there is a significant ( > 40%) 
component of CDM [41]. Since the neutralino is a CDM can- 
didate in our model in the case where there is an unbroken 
R-parity, mixed dark matter is a possibility. This does not 
diminish the advantage of a dark matter candidate which is 
necessarily warm. 

Discussion and Conclusions: We have shown that the Z2- 
singlino can account for the observed phase-space density of 
dwarf spheroidal galaxies, which may be evidence of non- 
singular cores. In the case where the % mass comes entirely 
from the Higgs expectation value, Mj is fixed by the % relic 
density. The observed abundance of dark matter implies that 
Mj£ is in the range 0.3-4 keV, which coincides exactly with 
the range required for % to act as WDM. Therefore dark mat- 
ter is necessarily warm in this model. The model can directly 
account for the phase-space density of dwarf spheroidal galax- 
ies when Tr « 10 — 100 TeV, while dynamical suppression of 
the primordial phase-space density allows smaller values of 
Tr to be consistent with dSphs. The free-streaming length is 
in the range 0.3-4 Mpc, the lower end of which may reduce 
the overproduction of satellites and loss of angular momen- 
tum observed in CDM simulations of galaxy formation. The 
small mass of % can be understood in terms of a large mes- 
senger mass, Ms ~ 10 10 GeV. Such a heavy S field might be 
identified with the messenger sector of gauge mediated SUSY 
breaking models. We will return to this possibility in a future 
study IhIIi . 

Depending on Tr and sin2|3, a range of primor- 
dial phase-space densities can be generated, with 
Q w 10 5 - 10 10 eV/cm 3 /(km/s) 3 when sin2|3 > 0.1 
and 10 2 GeV < T R < 10 5 GeV. This may allow the 
wide range of observed values of Q, ranging from order 
10 6 eV/cm 3 /(km/s) 3 in dSphs to 10 4 eV/cm 3 /(km/s) 3 or 
less in normal spiral galaxies l44ll . to be understood, for 
example by having the values of Q in dSphs close to their 
primordial values and the values in normal spirals suppressed 
by non-adiabatic evolution during structure formation. It has 



' Similarly, constraints on the WDM particle mass from globular cluster- 
based observations of the phase space density of the Fornax dwarf galaxy 
14211 cannot be directly applied to the present model. 
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been suggested that too many dSphs may be generated when 
the mean primordial Q is equal to that in dSphs ifTll l45ll . 
In that case a possible solution might be to have a mean 
primordial Q much smaller than that observed in dSphs, 
with dSphs then forming from a high phase-space fraction of 
the Y p articles in the low momentum tail of the distribution 

MM- 

In addition, strong constraints on % dark matter may be ex- 
pected from Lyman-a observations. However, existing con- 
straints on WDM masses cannot be directly applied as the % 
momentum distribution function due to thermal Higgs anni- 
hilation differs from those considered in existing studies. We 
will return to these issues in a future study l43ll . 

If /^-parity is unbroken in the MSSM then the model can be 



extended to a mixed dark matter model, with the /^-stabilized 
MSSM LSP providing CDM in addition to the Z 2 -stabilized 
% WDM. In this case WDM should easily be compatible with 
Lyman-a constraints. 

Testing the model at colliders will be challenging due to the 
small effective coupling of % to the Higgs bosons, with the hjx 
coupling being of the order of v/M s w 1CT 8 . However, the 
simple form of the superpotential may allow the model to form 
part of a more complete model which could have distinctive 
collider signatures. 
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